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UPt3 exhibits a paradoxical behavior in its antiferromagnetism. Namely, the neutron scatter- 
ing experiments surely detect the antifferomagnetic ordering below the transition temperature 
Tn ~ 5-fT, while no trace of the transition is seen in the NMR experiments and in the specific 
heat at Tjv- For resolution of this paradox, we propose the "SDW state" which oscillates with a 
finite frequency larger than those of the NMR and smaller than those of the neutron scattering. 
For demonstrating such possibility we seek the divergent point of susceptibility, and offer a 
preliminary calculation of the frequency dependence of "magnetization" . 
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The heavy fermion compound UPt3 has paradoxical magnetic properties. Neutron scattering 
experiments show the onset of antiferromagnetic order at TV= 5K0> with the presence of ordered 
moment equal to 0.02/i£ at low temperatures. On the other hand, NMR experiments cannot detect 
any trace of the ordered moment and the specific heat has no signal of the transition at all.B These 
experimental facts are considered to be well established. These apparently contradictory results 
may be understood naturally if the "antiferromagnetism" detected by neutron scattering is not a 
true order but fluctuates with some finite frequencies smaller than ~ 10 11 /s, characteristic of the 
neutron scattering, but greater than ~ 10 7 /s, which are those of the NMR measurements. However 
the correlation length, £, is about 300A which is so long that it is nearly equal to the resolution limit 
of the neutron scattering experiments. In addition, it can be observed over a wide temperature 
range from the "transition temperature" 5K down to zero temperature,! this fact indicates that 
the spin fluctuation does not obey the ordinary scaling law, £ oc \(T — T c )/T c \~ u . We consider that 
this apparent paradox indicates a novel phenomenon and we propose here as its solution that the 
"SDW state" oscillates with a finite frequency. 
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First we explore the possibility that the divergence of staggered susceptibility occurs at a finite 
frequency rather than at zero frequency. If it is possible, it strongly suggests the existence of the 
"order" which may oscillate with some finite frequency. The "order" we mean here is a kind of spin 
fluctuation which can be regarded as an antiferromagnetic state if observed for a very short time 
interval like the neutron scattering investigate, in other words, we may call it a "dynamical order". 
We search the origin of such an exotic state for the nesting property of the Fermi surface which 
causes strong frequency dependence to the mode-mode coupling coefficient The nesting property 
of the Fermi surface of UPt3 is partly supported by the band structure calculations and the dHvA 
effect .0* We calculate the frequency dependence of "transition temperature" using the mode-mode 
coupling scheme with a mode coupling coefficient which has a peculiar frequency dependence. 

On the basis of the itinerant-localized duality modeljli^ the dynamical spin susceptibility x(Q + 
q, \0J m ) for the coherent spin fluctuations is given as, x(Q + Qi ltjL, m) = Xs(Q + q, i^m)+X/(Q + ( /i luJ m)j 
where Xs{Q + (M^m) is due to the localized components of spins, Xf(Q + Qj '^m) to the itinerant 
fermions, and Q is the commensurate antiferromagnetic wave vector. 

The localized spin part Xs(Q + q,iw m ) is obtained using the mode-mode coupling theory as 

X7 1 (Q + Q,'u>>m) = Xso(Q + <l' iuJ m ) -U(Q + q,iuj m ), (1) 

where the dynamical susceptibility for the localized spin Xso(Q + Qi^m) without the coupling to 
the itinerant fermion is given as x^o* (Q + <7j ' 1UJ m) = Xo 1 fam) ~ J{Q + 0.) > where xo (i^m) is the local 
susceptibility and J(Q + q) the RKKY-type interaction among the localized component of spins. 

The fermion polarization II is given as II = IIo + IT where Ho is the bare electron-hole polariza- 
tion part and the IT is the mode-mode coupling term, which are given as 



IIo(Q + q,iu m ) = -g 2 T^G(p,ie n )G{p + Q + q,ie n + iu m ), 



n.p 

Z 2 



= -g 2 T V -rj; ; r (2a) 

t£ (ie n - £p)(ien + W m - Zp+Q+q) 

Ux(Q + q,iu m ) = -g 4 T ^ G aX s(Q + </, M m ,), (2b) 

m',q' 

where G(p,ie n ) is the Green function for the itinerant fermion, z is the renormalization amplitude, 
and g is the coupling constant between localized spin fluctuations and fermions and is related to 
the fermion coherent temperature zg ~ To® The mode-mode coupling coefficient (G4) is given as 
G 4 = 6G| - Gl where 

G l = T E E [e ^)G(p + Q + q,ie n + iu m ) (3a) 

n p 

xG(p + Q + q',ie n + iO m /) 5 
G l = T E E G ^ [ ^)G(P + Q + q;ie n + iLO m ) (3b) 

n p 
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Fig. 1. The Feynmann diagram for the mode-mode coupling coefficient, G4. The solid line represents the fermion 
propagator and the wavy line the spin fluctuation propagator. The dot represents the spin-fermion coupling 
coefficient. 



xG(p + q - q',ie n + iU m >)G(p + q' - Q,ie n - itl m ')- 

The polarization IIo, given by eq.(2a), is expanded for small \q\ and small u> as, Ho(Q + q,u>) ~ 
A r i ^ 1 [l — Afq 2 + iCu], where we assume H(Q,0) to be Np l with Np being the effective density of 
states. The coefficients Af and C are given as!) Af = afV F /[max(T, h)] 2 and C = Cf/max(T,h), 
respectively, where ay and cy are numerical constants, vp is the renormalized Fermi velocity, and h 
is the parameter measuring the deviation from the perfect nesting condition. The nesting property 
of the systems is reflected in the strong temperature dependence of the coefficient Af and C in 
contrast to the systems without nesting. So the retarded counterpart of Xs{Q + Q, i^m) is expressed 
as 

Nt 

Xs{Q + q,uj) = 2 1 2 — — , (4) 

A[k{lo) z + q z \ — \Cuj 

where k is the inverse of the correlation length, A = Af + a s Ja 2 N f where a s are numerical constants 
due to the localized spins, J is the energy scale of the RKKY-type interaction and a is the lattice 
const ant 

Due to the nesting property, G4 in eq.(2b) not only has a temperature dependence but also has 
a strong frequency dependence. Here we calculate the mode-mode coupling coefficient assuming 
that there exists a 'dominant frequency' iuj m as well as a wave vector Q. The mode-mode coupling 
coefficient, which is shown in Fig. 1, is obtained after performing frequency summations as, 

G 4 = z 4 g 4 y (th% - th%g ) L- (5) 



4T ^ K 2T 2T ' (iuj m + e p - t P+Q ) 2 

When the perfect nesting condition £ p +q = — £ p is fulfilled, we can see that G4, given by eq.(6), 
has asymptotic dependence as \oj m \ — > 00. In Fig. 2 we also show the calculated result of 
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frequency dependence of (6) with the model dispersion £ p = — 2t(cosk x a + cosk y a) - /ion the square 
lattice for various fillings. We can see the strong uj m dependence which drastically decreases the 
mode- mode coupling coefficient with asymptotically proportional to l/u^- This marked decrease 
in mode-mode coupling coefficient reflects the nesting property and is important for our discussion. 
Combining the asymptotic dependence of temperature we can take the interpolating dependence 



Fig. 2. Frequency dependence of the mode-mode coupling coefficient with the model dispersion £ p = — 2t(cos k x + 
cosfc H ) for various filling n. The bold line represents the slope of uj~ 2 for a guide to an eye. 



of G\ on T, h and u m , for simplicity, 



r 2 . 



where r = \/T 2 + h 2 and gj is a numerical constant from fermion contribution and is given as gj = 
35£(3) /8ir 2 for the 3D spherical band and F is the parameter which is introduced so as to reproduce 
the correct asymptotic form. We choose this interpolating form so as to reproduce the correct 
asymptotic form concerning dependence on uj m and T, and not to lead to unrealistic singularities 
when Cr4(iu; m ) is analytically continued onto the real axis, i.e., itu m — > u + \8. Considering that 
x(Q,iuj m ) takes real values and is as even function of u m , we may be able to use eq.(|6|) as an 
appropriate form. With this frequency dependent mode- mode coupling coefficient, we follow the 
conventional mode-mode coupling theory The self-consistent equation is derived from (Q) (2a) 
and (2b) as 

y(ko m ) = ^-( x7o \Q,Q) -n (Q,0)) 

+^2Gi{\uj m )^T^Xs{Q + q,iuJ m >) (7) 

A 1r n m , 
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Then we arrange it as follows, 



r 2 „ , , r 5 



y(iu) = 2/0 ., , „_ 2 + Fg{u) - ■ _ (8) 



1 + qt 2 (1+ar 2 



/ <ixx 2 / 
Jo Jo 



00 , 1 CJ 



#>/t _ 1 + x 2]2 + [_ 



/3(l+aT 2 )J 



where 

and i, r, h, and u; m have been normalized by dividing by iV^ 1 ~ ep, and i* 1 = gfNpz^g^q^^ir^cA 
is the constant which is proportional to the strength of mode-mode coupling. yo = 
(NF/af)[xsO~ 1 {Q, 0) — no(Q, 0)] is the measure of the distance from the magnetic phase boundary 
at T = 0. a = Aa s J/7rafeF is the normalized RKKY-type interaction and /? = atjlixcf is the 
numerical factor. For the 3D spherical band, a and (5 are given as a = (96-7r/7£(3)) J/e^ and 
/3 = 7C(3)/87r 4 , respectively, is the physical cut off wave number and is given as qs 3 = 6tt 2 /vo 
where vq is the volume per magnetic atom. y(u>) in the integral (8) is analytically continued from 
yiiLOm) onto the real axis in complex to plane. The self-consistent equation for y{iuS) is written as 

r 2 r 5 

y{iuj m ) = y — 2 + F 9{um)j7- 2\3 S ( T )> ( 10 ) 

1 + ar z (1 + aT z ) A 



where 



S(T) = (11) 
dxx 2 f duj- 



lo Jo e^-ny^+x^ + i^^r 

First we have to determine the parameter H(T) self-consistently from eqs.(|lO|) and (11). Then we 
determine the "transition temperature" Tjv(iu;) from the equation, XsiQ,^) 1 = 0, i.e. y(iu>) = 0, 
for finite frequency \ui for which Tjv(iw) is expected to be higher than Ttv(0). Namely, Tjv(iu;) is 
determined by the condition 

1 r 4 
y r + — u + Fg(u) E(T) = 0, (12) 

2irp (1 + a.T z y 

where the frequency \u) is analytically continued on the imaginary axis from the Matsubara fre- 
quency points iw m 's. 

The numerical solution of eq.(|i~2|) is shown in Fig. 3. We can see that the highest "transition 
temperature" Tjv is obtained at the finite frequency iluq. The abrupt disappearance of the transi- 
tion temperature in the high frequency region is due to the w-dependent structure of eq. (|l2|) , which 
may not be valid in this region. A realistic physical situation may be that the "transition temper- 
ature" once reaches its maximum at some frequency, iivq, and then gradually approaches zero with 
increasing frequency \tu. Hereafter we write the highest "transition temperature" as Tjv(iwo)- 
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The reason for such frequency dependence of TV is easily understood by inspecting the structure 
of eq.(|l2|). The mode-mode coupling term of eq.(fL^), which plays the role of reducing the transition 
temperature T^r, is suppressed by the factor g(u>), given by eq.(^), with increasing frequency uj so 
that the effect of increasing Tn overwhelms the damping term arising from the electron- hole pair 
excitations, the second term of eq.(p^), in some region of frequency. 

In Fig. 3 we also show the result with large damping (small value of (3 in eq.(p"2|) ), in which case 



Fig. 3. "Transition temperature" Tjv(itj) as a solution of eq.(|l2|). 



the transition temperature takes the highest value at zero frequency like the ordinary SDW case. 
Thus we obtain a possible structure of the magnetic susceptibility that has divergence at some 
finite frequency under certain conditions. We recognize this as an instability due to the emergence 
of the dynamical component of the SDW moment. So the state we consider below is different from 
the usual paramagnetic phase and in this context it can be said that it is a phase transition but 
the order in this state cannot be defined for a long time interval. 

By a simple argument we can expect that the "SDW moment" Mq has its peak at some finite 
value of real frequency uj. If the mode-mode coupling coefficient given by eq.(||), which corresponds 
to the fourth-order term in GL expansion, has its peak at uj = 0, the existence of zero-frequency 
component of "order parameter" involved more energy than the other components of frequency. In 
order to confirm the above prediction we calculate the SDW "moment" on the real w-axis. 

For simplicity, we retain only one frequency component of the "order parameter", Mq(w), and 
do not take into account the coupling between different frequency components of Mg(w) as the 
first approximation. In the presence of Mg(w), the self-consistent mode- mode coupling equation is 
given as 

r 2 F r 5 

y± = yor— t — 2+ t 9 ^tt, — 2\3( 3H ^ + 2H ll) 

1 + ar z 5 (1 + a.T z ) A 11 



6 



+ j-^^jSMMq^) 2 (13a) 

n = ^ + M ^M^ry <13b) 

where the subscripts _L and || denote the transverse and the longitudinal components of suscepti- 
bilities, respectively, and y± = NF/Aqf 3 x±(Qi UJ ). We should note that the damping term, linear 
in iu> vanishes when the wave vector of x±{ ( li UJ ) is just the antiferromagnetic wave vector Q in the 
presence of Mq(w). The E u (T),i/ =_L or || is given as, 

S*(T) = (14) 
-1 roo i 



f f c 
/ dxx 2 
Jo Jo 



duo- 



The magnetization Mq(lj) is determined from the equation, x±(Q, w > Mq(c<j)) _1 = 0, where the 
subscript _L denotes the transverse component of susceptibility, and this condition is rearranged 

as, 

y±(Q,u) = yoT+ — 2 uj 2 (15) 

+ f 9{uj) (TTqW (3H± + 2H||) + T9{uj)Mq{uj)2 = °- 

Here we have inserted the w 2 -term in eq.(15) which is the leading term in to of the susceptibility 
for q = Q and the coefficient D should be determined by the band structure. We can determine 
Mq(w) by solving the coupled equations (13) ~ (15). For simplicity we approximate Ey by Sj_. 
This treatment is valid at low values of Mq(w). Since we are detecting the frequency dependence 
of the "order parameter" Mq(w) near the "transition temperature", this approximation may be 
sufficient for our purpose. The plot of Mq(lj) is shown in Fig. 4. We can see that Mq(w) has large 
values at finite frequency than at lo=0. Furthermore, Mq(w) vanishes for low values of uj (the case 
t=0.1 in Fig. 4). Such a situation is expected. There exists no Mq(0), the static component, so 
that it is not a true SDW state. We should note that a finite value of Mq(w) is only meaningful 
in the time interval ~ 2tt/uj, so that the existence of Mq(lo) indicates that there exists a spin 
fluctuation mode whose "correlation length" is extremely long and it as if there was a true SDW 
order when observed in the time interval ~ 2tt/lo. The specific heat will not be changed by the 
appearance of Mq(w) unless Mq(0) appears because the specific heat jump in the usual phase 
transition mainly comes from the static component of magnetization. So far, we have assumed 
many parameters with arbitrary values. But we believe that the real behavior of UPt3 is just 
in the parameter region where this novel phenomenon can occur. The Plot of Mq(uj) has been 
obtained appears by assuming that there exists no coupling between the "order parameter" with 
different frequencies. The SDW "order" first develops at some finite frequency. With decreasing 
temperature, the other frequency component, including the static (zero-frequency) component of 
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Fig. 4. Magnitude of the "ordered moment" Mq(lu) on real u>-axis. 



the SDW state, should gradually develop in general. Indeed, the ordered component Mq(w) can 
induce other components through the coupling term such as Mq(w) 3 • Mq(— 3a;), where Mq(— 3w) 
is the admixture frequency component. However, as long as the coupling term maintains its form 
similar to eq.(||), the induced order Mq(— 3a; ) whose frequency is located in the frequency region 
where Tn(31o;) = 0, results in a considerable energy cost and hardly grows. When the temperature 
is decreased and the 'ordered' state develops sufficiently, the mode-mode coupling term no longer 
has the frequency dependence like eq.(^). Then the dominant frequency of the SDW moment may 
change from the value near the "transition temperature". Therefore, there is a possibility that 
the specific heat jump at T=20mKH' may be interpreted as that at this temperature the dominant 
frequency of the SDW moment has changed to the static component (i.e., o;=0) and it reveals its 
own state as the ordinary SDW state.B* 

In conclusion, existence of an "oscillating SDW" state can possibly be verified using a probe 
with time scale intermediate between neutron scattering and NMR: say, /^SR measurements. % A 
possibility of the "oscillating SDW" state can be seen in a recent numerical simulations for an 
extended Hubbard model with pair hopping terms© In the ordered state we can determined only 
the qualitative level near the 'transition temperature' and further investigation is required to say 
details of this new state of magnetism. 
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